arXiv:hep-th/9507156v2 28 Jul 1995 


October 1994 


LBL-37560 

UCB-PTH-95/28 

hep-th/9507T^ 


Convergent sequences of perturbative approximations for 

the anharmonic oscillator 
II. Compact time approach]] 

B. Bellet 
P. Garcia 

Laboratoire de Physique Mathematique 
Universite de Montpellier II-CNRS 
34095 Montpellier cedex 05 
and 

A. Neveu| 

Theoretical Physics Group 
Lawrence Berkeley Laboratory 
University of California 
Berkeley, California 94720 


Abstract 

We present an alternative pathway in the application of the variation 
improvement of ordinary perturbation theory exposed in [1] which can 
preserve the internal symmetries of a model by means of a time compact- 
ification. 


*This work was supported in part by the Director, Office of Energy Research, Office of 
High Energy and Nuclear Physics, Division of High Energy Physics of the U.S. Department of 
Energy under Contract DE-AC03-76SF00098 and in part by the National Science Foundation 
under grant PHY-90-21139. 

lneveu@lpmsun2.lpm.univ-montp2.fr. On sabbatical leave after Sept. 1®^ 1994 from Lab¬ 
oratoire de Physique Mathematique, Universite de Montpellier II-CNRS, 34095 Montpellier 
cedex 05 






Disclaimer 


This document was prepared as an account of work sponsored by the United States Gov¬ 
ernment. While this document is believed to contain correct information, neither the United 
States Government nor any agency thereof, nor The Regents of the University of Galifornia, 
nor any of their employees, makes any warranty, express or implied, or assumes any legal 
liability or responsibility for the accuracy, completeness, or usefulness of any information, ap¬ 
paratus, product, or process disclosed, or represents that its use would not infringe privately 
owned rights. Reference herein to any specific commercial products process, or service by its 
trade name, trademark, manufacturer, or otherwise, does not necessarily constitute or imply 
its endorsement, recommendation, or favoring by the United States Government or any agency 
thereof, or The Regents of the University of California. The views and opinions of authors 
expressed herein do not necessarily state or refiect those of the United States Government or 
any agency thereof, or The Regents of the University of California. 


Lawrence Berkeley Laboratory is an equal opportunity employer. 



1 Introduction. 


In the preceding paper we presented a study of large orders of a variational 
method which provides convergent sequences of approximations for the quantum 
anharmonic oscillator. In this method, we introduced a mass term (with a mass 
a;) in the free Lagrangian and subtracted it in the interaction part. Then, we op¬ 
timized at any order in terms of this variational parameter lo to get approximate 
values of different quantities. 

We have in perspective the use of such a method in quantum held theory. 
For fermions like those of the Gross-Neveu model ||^, or QCD quarks, intro¬ 
duction of a mass term causes a priori no problem, but for gluons it would 
break gauge invariance, and there is no obvious modihcation of the momentum 
dependent part of the gluon propagator and the structure of the ultraviolet in- 
hnities of the theory. However, a space compactihcation might provide some 
good approximation of quark conhnement, like in the M.I.T. bag model 

In the M.I.T. bag model, the bag is a fundamental ingredient of the theory, 
which hence is not QCD, while we would advocate here an approach in which 
it would be added to the free gluon and quark action in the form of a modifi¬ 
cation of their propagators, the subtraction of this modification being done in 
perturbation theory. While formally when going to all orders one would have 
done nothing, and still be dealing with QCD, a hnite given order would have to 
be optimized in the size of the bag as the closest one could get to the inhnite 
order case. The advantage would be that different physics would be involved, 
hopefully closer to that of the full theory, with conhnement, etc... 

Compactihcation as an intermediate step in calculations in infrared diver¬ 
gent held theories has already been considered by several authors , with 
reasonable success. These authors usually treat the theory in a spatial box 
of length /, and extrapolate the results to / — oo. We are here advocating 
a diherent, if related, strategy, in which the box is in some sense subtracted 
perturbatively, and the calculation at a given order in perturbation theory op¬ 
timized with respect to the size of the box, this being supposed to mimick as 
well as it can the true all-orders theory, where formally there is no dependence 
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at all on the size of the box. 


In this paper, we try ont this idea on the case of the anharmonic oscillator, 
choosing an action with compactihed time as our starting point. We hnd a 
conhguration for the optima and the approximate values of the ground state 
energy very similar to those of the harmonic approach of this problem [Q. In 
the second section, we study the interaction part of the action to be used in order 
to formally leave the total action of a pure anharmonic oscillator unchanged. We 
make a choice which allows relatively simple perturbative calculations starting 
from a theory governed by a compactihed action. The perturbation theory for 
such a compactihed action is considered in the third section. In the fourth 
section, we establish the variational-perturbative expansion of the mean value 
of the Hamiltonian and optimize it in terms of our variational parameter (the 
size of the box) up to order 16. We interpret the conhguration of the extrema 
using the large order behavior of the expansion, which is qualitatively much the 
same as in the harmonic approach of the previous paper, to which we compare 
our results. They turn out to be less accurate numerically for similar orders in 
the approximation, which rehects the fact that approximating an anharmonic 
potential by a compact time dimension is much worse than approximating it by 
a harmonic potential. Nevertheless, the convergence properties are qualitatively 
very similar. 


2 Time compactification of a pure anharmonic 
oscillator. 


The Euclidean action for a pure anharmonic oscillator writes: 


S 





dt . 


( 1 ) 


Severe infrared divergences prevent any perturbation in powers of A. A way 
to tame these divergences is of course to introduce an term as was done 
in the previous paper. Here, we would instead like to replace in some sense the 
free term of the action of equation (1) by the compactihed free action S^{t): 
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( 2 ) 


= /l ^ dt , 

and rewrite the total action S as: 

S = S^{t) + S\t). (3) 

Then, we will be able to perform a perturbative expansion of any quantity in 
powers of the interaction term S\t) around the free term S^{t). To do so, let 
us consider the Euclidean action St of a time compactified pure anharmonic 
oscillator 


IL'T' 
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^ + ^0^ ] dt. 


( 4 ) 


With St, one can compute in perturbation theory in powers of A at any 
given T (provided one uses appropriate boundary conditions, see below). To see 
how S and St may be related, we perform the following change of variable: 


t 

a; = Ttan— ; dx = {I S—) dt, 
1 1 ^ 


under which the action St becomes: 


/ OO T \ T 

jd + (4^) + 4 77-^^“ I 


(1 + fj) 

On the other hand, we can rewrite the action S as 

1 


( 5 ) 


( 6 ) 


^ = 


roo r 1 o A 

I - (1 + ^ (1 - e)) (d,(Pf + - -^- 

-ooL2^ ^ 4:1 + 4 (l-e) 


(4 dx, (7) 


where, strictly speaking, 5 = 1. Comparing (7) with the expression of the 
compactihed action (6), we see how to relate S and St'- following the procedure 
of the previous paper, one computes the quantity of interest to some order in 
A using St, sets T = r(l — 5 )“^/^, expands in powers of e to the same total 
order as in A, sets 5 = 1 and optimizes the result with respect to r. This should 
provide us with an approximate value of the same quantity computed with S. 


In more detail, we can consider the action (7) as: 


S = S\t) + S\t) 
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with S^{t) given by 


and 



S\t) = 


1 X 

'o ^ - 
Z T 


{9x4')^ + 


4 1 + fi (1-e) 


dx. 


( 8 ) 

(9) 


The order in perturbation will contain all the terms proportional to 
n + m < p. Note that in this approach, we do not take into account 
the complete quartic term of (4) but that we reconstruct it more and more 
accurately as the order of perturbation increases. This feature does not seem 
to us to carry any pathology, but simply follows from the choice of St as our 
compactified action and the physics it contains. 

This new choice for a free kinetic term parameterized by r in S' will give rise 
for every hnite order p of perturbation to an explicit dependence on r. Because 
the total action S remains unchanged at e = 1, one expects this parameter r 
to become irrelevant at infinite order of perturbation. As already explained, 
this suggests to fix this parameter by looking for the regions where the result to 
order p is stationary in r. 

Let us now discuss briefly the boundary conditions in the action of equation 
(4). Periodic conditions (()-t = 4t sxe unsuitable for perturbation theory, as 
the Gaussian kernel has an isolated zero eigenvalue. Similarly for Neumann 
boundary conditions. Hence, we shall use Dirichlet boundary conditions in which 


4±t = 0 . ( 10 ) 

which preserves the discrete 0 —0 symmetry. With this choice, all infrared 

problems disappear at hnite T. 

In order to compare the results obtained by this compact time method 
to those of the harmonic one [Q, we would like to consider the ground state 
energy. However, the method we described above implies calculation using a 
compactihed action which is therefore not invariant under translation in time. 
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Then, the ground state energy is not a meaningful concept and we choose instead 
to consider the mean value of the Hamiltonian operator at some time to'. 

< K > ( 11 ) 

The variation over the size of the box T provides some new approach of the limit 
T —cx) (substituting r —T(1 — expanding in powers of e and taking 

e = 1). In the limit of infinite order of perturbation, the box should disappear 
and we should recover the value of the ground state energy when calculating 

< 4o >• 


3 Perturbative expansion for the compactified action. 


In the Euclidean path integral representation, the mean value of the Hamiltonian 
operator at time to in the box writes: 


where 


< Ht, >= 


< 0; T I Ft J 0; -T > _ 
<0;T|0;-T> ’ 


( 12 ) 


< 0;T I 0; -T > = / D4> Ht, , 

J 0±T—0 

< 0; T I 0; -T > = [ , 

J 0±T—0 

-St[ 4>] = - +^(j)t]dt , 

and Ht, = ■ 


3.1 The partition function and the mean value of the 
Hamiltonian. 

Dividing by the partition function Z{T) in (12) cancels all the disconnected 
Feynman diagrams. We rewrite it: 

Z{T) =< 0; T I 0; -T >=< 0 | f/sT | 0 > 
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Here, in the Schrodinger representation, the evolution operator Ut of a quantum 
system described by H appears. It satishes the following equations: 

-jUt = HUt and Uo = 1. (13) 

Then, taking the derivative of the partition function with respect to the size of 
the box, we get: 


d , , , d - , 

^Z(T) = < 0 |^£/..| 0 > 

= —2 < 0 I HJJ 2 T I 0 > 

= -2 < 0;T I .Hi I 0;-T >, 

so that 

Tz(r) = -2 < + jit e [-t-.t], (i4) 

Thus, the mean value of the Hamiltonian operator does not depend of the time 

when it is taken in the box. Moreover, its perturbative expansion in powers 

of T is obtained from the one of Z{T). However, the variational procedure 

implies that at the order p in perturbation, the order contribution which is 

proportional to A”'r'^(l — has to be expanded in powers of £ up to order 

p — n. This power of A does not take into account the A present in front of 

the quartic part of the Hamiltonian. But the two sides of equation (14) identify 

themselves in perturbation, mixing that A with the perturbative one So, we 

have to know kinetic and potential contributions for each order in the expansion 

of dZ{T)/dT. Nevertheless, we will use this equality to compute at any order 
" 2 

the value of < {dtcfy^ > in term of < 0^ > and Z{T). We decided to compute 
the mean value of the Hamiltonian where its kinetic and potential contributions 
do not much vary and are the least dependent on the size of the box, i.e. in the 
middle of the box at t = 0. 


*For example, in the perturbative expansion of dZ{T)/dT, the zeroth order contribution 

» 2 

of < bt > is mixed with the first order of < > ■ 
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3.2 The perturbative calculation at any order. 


In order to develop a feeling on how the variational method works, we have to 
reach high orders of perturbation, i.e. calculate a large number of coefficients. 
The Feynman diagrams approach implies the calculation of an increasing num¬ 
ber of graphs. Instead, we prefer a perturbative expansion of the evolution 
operator closer to the usual Rayleigh-Ritz method. Expanding Ut in powers of 
the coupling constant A, we get the following expansions for Z{T), < 0q > and 

< (<9*0)0 >: 

1 00 

<0;T|0;-T> = ^ (-A/4)"T3- Z„, 

V n=0 

1 00 

<0;r|.^J|0;-r> = ^ ^ (-A/4)”r»+" Q„, 

V -L n=0 

'■'2 1 

<0;T\{a,4>)„\0-,-T> = ^ ^ (-A/4)”r»-‘ P„, 

V 4 n=0 

with 



The zero subscripts appended to ket vectors mean that the corresponding expec¬ 
tation values are computed using free propagators with the boundary conditions 
discussed in section 2. 


In the above expressions, we have extracted the T dependences of the co¬ 
efficients Zn, Qn and Pn. In appendix, a systematic calculation of the Z„’s and 
Qn’s is performed with the help of a recursive formula. The values of the Pn’s 
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are deduced from these using (14), which gives: 


Pn — (3u — 1/2)Z„ — 2Q„_i. 


After performing the division by Z{T) order by order in perturbation, the 
expansion of the mean value of the Hamiltonian writes: 

< H„ >“= D-5(x>” P" 

where the c superscripts of P and Q refer to the contribution of the connected 
diagrams only. We have computed these coefficients up to order 16, and we give 
here the first few of them: 


QS 

Q5 

Qi 

Qs 

Ql 

Qs 


3 

4 

10 

15661 

700 

3798833 

19250 

634428694707 

297797500 

10113181264708 

372246875 


P5 

Pi 

pc 

■*^^3 

pc 

pc 


-1 

2 

10 

349 

175 

2363729 

250250 

14511295339 

223348125 
8165883862419 
14145381250 ' 


(16) 


4 The results and their interpretation. 


Starting from expansion (15) at order p, we perform the replacement > 
t‘'{1 — Expanding in powers of £ up to order p — n, we get the required 

expansion, which then has to be optimized with respect to our variational pa¬ 
rameter r: 


= '^1—P'^ r(p -I- n/2 -I- 1/2) 

° ^ ^ 2 ^ 4 ^ " r(3n/2 + l/2)r(p-n + l) 

, QC r(p -F n/2 -F 2) 

+ 4 I 4 J r(3n/2 + 2)r(p-n + l) 


(17) 


The values of this expression at its optima in r will serve as variational estimates 
of the value found in [^: 


Eexact ^ yi/3o.420805 ... 


(18) 
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Such a polynomial (17) has a number of optima which increases with the order 
of perturbation p. Most of them are complex and give a small imaginary part to 
the estimated value of the energy as in the harmonic approach [^. All of these 
estimated values exhibit the expected factor. We take A = 1 and compute 
all the optima up to order 16. Results and comments follow in the next section. 

4.1 Solutions of our variational problem. 

(r) is plotted in hgures |I| and versus the variational parameter for several 
even and odd orders. 

Since r is supposed to be irrelevant at inhnite order of perturbation, as 
expected, as p increases, the curves (r) flatten around the exact value over 
an increasingly large range of values of r as p increases. 

We can follow on these hgures the evolution of the optima with the orders 
of perturbation. The hrst minimum to appear provides a rather poor estimate 
{Eq ~ 0.34...) but moves to the left of the hgure as the order increases, and 
is replaced by a much more accurate maximum [Eq ~ 0.43...). It is clear 
that a third extremum, a minimum, appears with a still more accurate value 
Eq ~ 0.417..., and, from the trend of the odd orders, another, presumably even 
better maximum would appear if one would push the calculation to a couple 
more orders. 

The polynomial Eq^^ to be optimized has an increasing number of station¬ 
ary points as the order of perturbation increases. Plotting real optima versus 
the order of perturbation up to p = 16 (hgure 3) reveals how they arrange them¬ 
selves: Optimal values of r belong to families characterized by the straight lines 
in hgure 3 with slopes decreasing proportionally to the inverse of the square root 
of the order p. Such a family is created at hrst order, followed by another one at 
fourth order. One reasonably expects other families to appear at orders greater 
than 16, providing some even more accurate values. Figure 4 plots the location 
in l/r^ of the all the optima up to order 16 in the hrst quadrant of the complex 
plane. One thus easily follows each family up to order 16. 
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Such a regular behavior of all these optima is explained in the next subsec¬ 
tion by considering the asymptotic expression of Eq^'^ (r) for large p. 


4.2 The large order behavior. 


For large p, using the Stirling formula, we hnd: 






n=0 


r(3n/2 + 1/2) 


(^pl/2)3n-l ^ 


4' 4 




r(3n/2 + 2) 


{TP 


l/2\3n+2] 


(19) 


where we have taken A = 1. So, at large orders, the function to be optimized 
becomes a series in powers of We voluntarily omitted its upper bound 

in n. Using onr knowledge of the first coefficients, we establish empirically the 
following asymptotic behaviors: 

PC ~ 4 2 ”r(n), 

Q“ ~ I 2 "r(n + 2 ). 

With these behaviors, the term of the series decreases like and the 

series has an inhnite radius of convergence. For large orders p, the search for 
stationary values of E^\t) in r is nothing but the search for those of the 
series (19) in U = l/(r^p). Then, real extremal values of Y are the slopes 
of the families on hgure 3. Using the optima already calculated, we can then 
extrapolate the r optima and the corresponding valnes of the energy for large 
orders. 


The optima of the real family appearing at first order behave asymptotically 
as: 

l/(^opi)2 _ + 2.31156 p + . (20) 

p 

The corresponding valne for the energy tends toward a rather bad approximation 
from above: 

eT ^ 0.335 + . (21) 

p 
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( 22 ) 


For the second real family which appears at the fourth order 

l/(^opt)2 _ Q_4g5^ ^ 0.3578 p + , 

P 

the corresponding energy tends towards a better value from below: 

~ 0.434 - . (23) 

P 

The asymptotic value provided by a family is more accurate when the family 
appears later in perturbation. We also verify that the slopes tend to decrease. 
This corresponds to the successive extrema of (22) appearing at smaller and 
smaller values of Y, corresponding to those of the function Eq^{Y) 


°° 1-1 

KiY) = EI-o(—)”P 

n=0 


2^4' ”r(3n/2 + l/2) 


4^ 4 ' 




r(3n/2 + 2) 


y-(3n-l)/2 _|_ 


-^-(3n+2)/2j 


(24) 


Using the available coefficients up to order 16 is sufficient to see that this function 
has a minimum: 


Y ~ 2.2629 ; - 0.3315. 

One recognizes the slope and the asymptotic energy of the hrst real family within 
a few percents. Using the empirical asymptotic behaviors of the coefficients P 
and Q would extend the range of Y where one could compute Eq^, revealing the 
asymptotic properties of the other families. 


5 Conclusion 

The results presented in this paper are qualitatively very similar to those ob¬ 
tained in the harmonic approach [^: In the same way as in the harmonic ap¬ 
proach, the optima arrange themselves in families which are understood using 
a large order behavior. The set of these families provide a sequence of approxi¬ 
mate values which converges to the exact one. However, for a given perturbative 
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order, our estimated values were more accurate in the harmonic approach. In¬ 
deed, for the hrst family, we obtain an estimated value of the ground state 
energy with a precision of 2.10“^ in the harmonic case and of 2.10“^ in the 
present approach.This is presumably due to several causes: in the compact time 
case, there is no time translation invariance, hence no Hamiltonian, and we only 
compute the expectation valne of an operator which becomes the Hamiltonian 
in the large time box limit. Furthermore, as explained in section 2, our per¬ 
turbative expansions do not take into account the complete 0^ interaction term 
but only reconstruct it more and more accnrately as the perturbative order in¬ 
creases. Finally, a compact Enclidean time is clearly quite far from the physics 
of an anharmonic oscillator, much farther than a harmonic approximation. One 
may even consider it qnite remarkable that the procednre nevertheless seems to 
converge in the same manner, albeit much more slowly. One could thus con¬ 
template using snch a compactification, which could involve only the spatial 
coordinates or both space and time coordinates as a gange invariant variational 
approach to gauge theories. 


6 Appendix: Calculation of coefficients. 
6.1 Calculation of the Z„’s 


CX) _ \ 

<0;T|0;-T> = ^(— 

n=0 ^ 

Zn = f dti...dt„<O;l|0i .. .0i|O;-l>o 

For each time ti, ^ 2 , • • •, , one uses the closure relation in position space. 

The (j) operator is diagonal in this basis. One thus rewrites the integrand in the 
following way: 

< 0; -1 >0 = / dxn... dxi < 0; 1 | Xn, >o 

J —OO 

X < Xn]tn \ . . .\Xi]ti>o x\ < Xi] ti \ 0; -1 >o . 
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For t' > t, the free propagator writes: 

1 


< y; t'\x] t >= 


— t) 2(F t) 


{y - xy 


This product of expectation values of operators is a Gaussian function in terms 
of Xi, a; 2 ,..., a;„.One thus performs integrals over all positions simultaneously to 
obtain: 

< >0 = exp(^J^. A J) Ij^o 


2v^ 

^ (f. ,,,dy — _1 ^ 

(2n)!42' 


{r. A. jy ^, 


where 


Aij = (ti + 1)(1 - tj) for j > i {tj > ti), 

n 

detA = 2”"^ (fi+i - fi). 

i=0 

4„ is the fourth derivative with respect to Now, we build a recursive proce¬ 
dure which performs the four derivatives with respect to jn (4n) together with 
the integral over G [fn-i! !]• To do so, one will use the following identity: 






with 


Let us dehne: 


+p j=o 


r'p _ 

” p\{n — p)\’ 


fl{n,m,a,P)= f dti... f dt^ {J.A.jy {J.A) 


I J^O 


X (1 + tAr)“(l — ^at)^ I j=0 


= 4! 


\fdh... /' dtn-1 4 ... c\ (E 4 2-4 {JJUy 

“'iw-2 AT-I 

V-1 
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where: 


JJ^ = ji^ikjk and ^ ^ + 

]\f i,k=l JV *=1 


So that r2 satishes: 


4 4+a—i 

Vt{n,m, a ,/3) = 4lJ^Ai(n,m) 6ij(a,/3) 

i=0 j=0 

X n(n — i, m + 2i — 4,4 + a — i — j, i + /3 + 1) 


Mn,m) Ci—‘Cf2-‘ 


^i+/3 

_ '^i+a+P—j 


Q -(a 3) = _ 

CiP,, 5 + « + /3 


with the boundary conditions: 

r2(n, m, a, P) = 0 as soon as n or m is negative, 

r 2 ( 0 , 0 , a , P ) = ^ 11111^(1 + ^)“(1 ~ 2 ^, 

for the last integration to be performed between —1 and 1. The coefficients 
of the perturbative expansion are: 


6.2 Calculation of the Q^’s. 


<O;T|0^|O;-T> = 


1 \ 

d= E(-r)’V=’"+"Q, 

^ „.o 4 


E Q^-i 

p=l 


QLi = 




dt ^ • • • dtp • • • dt-yi 
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The preceding method must now be adapted to a slightly different domain of 
integration over the times (now tp = 0). We will consider four different regions 
in order to build the recursive dehnition of a new function flp{n, m, a, /5); 

rO pO pi pi 

/ dti... dtp-i “(tp = 0)” / dtp+i ... dtn 

J— 1 '^tn —2 V ^0 Jtn—1 

(4) ^ ^ (2) (1) 

These four different domains of integration can enter the preceding frame.The 
quantity (2n + m)/4 = n,..., 1 counts the remaining iterations to be performed 
and decreases by one unit at each step. One then builds the Op’s for each of 
these steps, writing: 


n' = n — i , 
a' = a — i — j , 


m' = m + 2i — 4, 

= i + /3 + 1. 


Sums over i and j are respectively performed between 0 and 4, and 0 and 4+q;— i. 


1. When (2n + m)/4 is greater than p + 1: Op is just as O: 

Qp{n,m,a,P) = ^ (3) Qp{n ,m',a', (3'). 

i j 

2. When (2n + m)/4 equals p + 1: Here tp, the lower bound of the integral 
is driven to zero: 

(l + tp)“ (1 — tp)^ \tp=o = 1 = (1 + tp)° (1 — tp)*^ 

Then, choosing a', (3' = 0, one can incorporate this case in the preceding frame: 

Vtp{n,m,a,(3) = Ai(n, m) 0ij(a,/3) Op(n', m', 0, 0) 

* i 


3. When (2n + m)/4 equals p: At this step, no integral has to be taken over 
tp so that there is no sum over j and the exponents a' and f3' remain zero: 

Op(n, m,0,0) = ^ Aj(n, m) Op(n', m', 0, 0) 
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4. When (2n + m)/4 is smaller than p: The domain of integration is 
obtained subtracting case (2) from case (1): 





dti ; 


Qp{n,m,a,p) = ^ ^ Ai(n,m) 0ij(a,/3) 

* i 

X {Qp{n , m', a', j3') — Qp{n',m',0,0)). 
This together with the boundary conditions 

flp{n, m, a, P) = 0 as soon as n or m is negative, 


fip(0,0,a,/3) = (5(a) 2^. 


dehnes Vtp for every p. 

The coefficients Qn-i are then: 


Qn—1 


1 1 1 
2v/i(2n)!4^ 


^ fip(2n, 0,0,0). 

p=i 


The values of the Pn’s are then deduced using formula (16) which writes in 
terms of the other coefficients: 

P, = (3n - 1/2) - 2 Q„_i Vn. (26) 
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Figure 1: E^\t) versus r for p =4, 6, 8 and 10. 



Figure 2: E^\t) versus r for p =7, 11 and 15. 
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Figure 4: All the optima in the complex plane up to order 16 





